Origin of Jahn- Teller distortion and orbital-order in LaMnOa 
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The origin of the cooperative Jahn- Teller distortion and orbital-order in LaMnOa is central to the 
physics of the manganites. The question is complicated by the simultaneous presence of tetragonal 
and GdFeOs-type distortions and the strong Hund's rule coupling between e g and t2 g electrons. To 
clarify the situation we calculate the transition temperature for the Kugel-Khomskii superexchange 
mechanism by using the local density approximation+dynamical mean-field method, and disentangle 
the effects of super-exchange from those of lattice distortions. We find that super-exchange alone 
would yield Tkk ~ 650 K. The tetragonal and GdFe03-type distortions, however, reduce Tkk to 
~ 550 K. Thus electron-phonon coupling is essential to explain the persistence of local Jahn- Teller 
distortions to > 1150 K and to reproduce the occupied orbital deduced from neutron scattering. 

PACS numbers: 71.10.Fd, 71.10.-w,71.27.+a,71.10.Hf, 71.30.+h 



The insulating perovskite LaMnOs is the parent com- 
pound of the colossal magneto-resistance manganites [1] 
and it is considered a text-book example of a coop- 
erative Jahn- Teller (JT) orbitally-ordered material [2j. 
Two distinct mechanism have been proposed to explain 
the cooperative distortion: many-body Kugel-Khomskii 
(KK) super- exchange (SE) and one-electron electron- 
phonon (EP) coupling [4]. Determining the relative 
strength of these mechanisms will provide a measure of 
the importance of strong correlation effects for the orbital 
physics in the manganites. Unfortunately the situation is 
complicated by the simultaneous presence of tetragonal 
and GdFe03-type distortions as well as a strong Hund's 
rule coupling between the Mn e g and tig electrons. 

In LaMn0 3 the Mn 3+ ions are in a t\ g e g configuration. 
Due to strong Hund's rule coupling the spin of the e g elec- 
tron is parallel to the spin of the ti g electrons on the same 
site. Above T/v = 140 K the spins on neighboring sites 
are disordered [5j . The crystal structure is orthorhombic 
(Fig. [lj. It can be understood by starting from an ideal 
cubic perovskite structure with axes x, y and z: First, a 
tetragonal distortion reduces the Mn-0 bond along z by 
2%. The La-0 and La-Mn covalencies induce a GdFcC>3- 
type distortion 0, 0] resulting in an orthorhombic lattice 
with axes a, b, and c, with the oxygen-octahedra tilted 
about b and rotated around c in alternating directions. 
Finally, the octahedra distort, with long (/) and short (s) 
bonds alternating along x and y, and repeating along z 
[l-[lll. This is measured by 8 JT = (I - s)/((l + s)/2). 
The degeneracy of the e g orbitals is lifted and the oc- 
cupied orbital, \9) — cos||3z 2 — 1) + sin||a; 2 — y 2 ), is 
~ |3Z 2 — 1), i.e. it points in the direction of the long 
axis. Thus orbital-order (OO) is d-type with the sign of 
9 alternating along x and y and repeating along z. At 
300 K the JT distortion is substantial, Sjt = 11%, and 
9 ~ 108° was estimated from neutron scattering data Q. 
Above Too ~ 750 K a strong reduction to Sjt = 2.4% was 
reported U UH , accompanied by a change in 9 to ~ 90° 




FIG. 1: (Color online) Structure of LaMnOa at 300 K 0]. The 
conventional cell is orthorhombic with axes a, b and c, and 
contains 4 formula units. The pseudo-cubic axes (left corner) 
are defined via a = (x — y)(l + a), b = (x + y)(l + P), and 
c = 2z(l + 7), with a,/3,7 small numbers. For sites 1 and 
3 the long (short) bond I (s) is ~ along y (x), vice versa for 
sites 2 and 4 (d-type pattern). All Mn sites are equivalent. 
The symmetries that transform them into a site of type 1 are: 
x <H- y (site 2) z — > — z (site 3), x y, z — > — z (site 4). 



Recently this was, however, identified as an order-to- 
Tfj,[llJ: Due to orientational disorder, 



disorder transition 1 



the crystal appears cubic on average, while, within nano- 
clustcrs, the MnO@ octahedra remain fully JT distorted 
up to T JT >1150 K [iH. 

Model calculations based on super-exchange alone can 
account for d-type order, but yield, for the classical 



ground state, 9 ~ 90° [13j|. Models of electron-phonon 
coupling in simple cubic perovskites instead give ~ 120° 
[J]. To explain the observed ~ 108°, one might thus con- 
clude that both mechanisms are of similar importance 
0- Such models are lacking, however, a realistic descrip- 
tion of the crystal and the calculated 9 is sensitive to the 
choice of parameters 0, Hi- LDA+J7 calculations yield 
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9 = 109° and show that Coulomb repulsion is fundamen- 
tal to stabilize the Jahn- Teller distortions in the ground 
state [l5j]. This might be taken as evidence that Kugel- 
Khomskii super-exchange is the dominant mechanism, 
and electr on-phon on coupling, enhanced by electron lo- 
calization 
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16], merely helps. On the other hand, re- 
cent semiclassical many-body calculations for model cu- 
bic perovskites indicate that electron-phonon coupling is 
essential to explain orbital ordering above 300 K [l7j . 

While it is not obvious how well LDA+J7 or semiclas- 
sical approaches capture the many-body nature of the 
KK super-exchange, it seems clear that the inclusion of 
the real crystal structure is crucial [3, H [HI 0|. The 
tetragonal and GdFe0 3 -type distortions result in a siz- 
able narrowing of the e g bands [1, 0, HO] , likely chang- 
ing the relative strength of super-exchange and electron- 
phonon coupling. Since, in the presence of a crystal- 
field. Coulomb repulsion suppresses orbital fluctuations 
jfiiHf . they may even compete with SE and EP coupling. 
To identify the driving mechanism for orbital-order in 
LaMnC>3, it is thus mandatory to account for both the 
realistic electronic structure and many-body effects. To 
understand the mechanism one has to disentangle the 
contribution of KK super-exchange from that of the JT 
or the GdFeCvtype and tetragonal distortions. 

In this Letter, we do this by calculating directly the 
Kugcl-Khomskii super-exchange transition temperature, 
Tkk, with and without tetragonal and GdFe03-type dis- 
tortions. We adopt the method used successfully for 
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FIG. 2: (Color online) Right: LDA+DMFT spectral function 
for the room temperature structure Rn for different U. ft (-11) 
indicates states with e 9 spins parallel (antiparallel) to St 2g ■ 
Left: k-resolved spectral function for U = 5 eV. 



orbital m! on site i' , obtained ab-initio by down-folding 
the LDA bands and constructing a localized e g Wannier 
basis. The on-site terms i — i' give the crystal-field split- 
ting. U and J are the direct and exchange screened on- 
site Coulomb interactio n |25j . We use the theoretical 
estimate J = 0.75 eV [26| and vary U between 4 and 
7 eV. The Hund's rule splitting was estimated ab-initio 



to 2JS t2g ~ 2.7 eV 0. We solve © using DMFT ^ or 



KCuF 3 [2l| , based on local-density approximation + dy- 
namical mean-field theory (LDA+DMFT) [Hj]. 

First, we calculate the electronic structure ab-initio us- 
ing the TV'^-order muffin-tin orbital method. Since the 
Hund's rule energy gain is larger that the e g -ti g crystal- 

and the e g bands 



bands are 



and 



field splitting, the ti g 

^-filled; the three t2 g electrons behave as a spin Sf 2 
couple to the e g electron via an effective magnetic field 
h = JSt 2g - In the paramagnetic phase (T > 1^=140 K) 
the ti g spins are spatially disordered. The minimal model 
to study the KK mechanism in LaMnOa is thus [24| 

H = V t 1 ' 1 ' u M ',c t c , , 

/ j m,m' <r,cr f ima v m' o ! 



cellular DMFT (CDMFT) and a quantum Monte Carlo 
[28| solver, working with the full self-energy matrix S mm ' 
in orbital space Q . The spectral matrix on the real-axis 
is obtained by analytic continuation 29]. 

We consider several structures: (i) the room tempera- 
ture structure, Rn with Sjt — 11%, and a series of hy- 
pothetical structures, R<s JT , with reduced JT distortion 
Sjt, (ii) the (average) structure found at 800 K, R|°4 K , 
which has a slightly larger volume than Rn and a smaller 
GdFe03-type distortion, and (hi) the ideal cubic struc- 
ture, Iq, with the same volume as Rn. For all structures 



we find that at each site the e„ spins align to St 



We 



calculate the orbital polarization p as a function of tem- 
perature [3(| by diagonalizing the DMFT (or CMDFT) 
occupation matrix and taking the difference between the 
occupation of the most (\8)) and least (\8 + it)) filled 
orbital. To test the t2 g spins picture we perform calcula- 
tions for the 5-band (e g + t% g ) Hubbard model (23|. We 
find that it holds even at high temperatures. 

For the 300 K structure (Rn) the band- widths are 
W t2g ~ 1.6 eV and W Cg ~ 3.0 eV. The e g states split 
by ~ 840 meV, in good agreement with experimental 
estimates [3l|. The lower crystal-field state at site 1 is 
|1> = 0.574|3z 2 - 1) + 0.818|x 2 - y 2 ). We find an insu- 
lating solution in the full range U = 4 — 7 eV (Fig. 
The Mott gap E g is ~ 0.6 eV for U = 4 eV, and in- 
creases almost linearly with increasing U. For U — 5 eV, 
suggested by recent estimates @, [32| , the Hubbard bands 
are at ~ —1.5 and 2 eV. In addition there is a broad fea- 
ture around 5 eV due to e g states with spin antiparallel 
is the LDA hopping integral from orbital m on site i to to the randomly oriented t2 g spins. These spectra are 
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c \mo creates an electron with spin a =ft,ft in a Wan- 



nier orbital 

J 



\x 2 - y 2 ) 



or |3z 2 — 1) at site i, and 
Uima = c\ ma c im(J . ft (ft) indicates the e g spin parallel 
(antiparallel) to the ti g spins (on that site). The ma- 
trix u {u^p, — 2/3 for i ^ i', u^ a , = 5 CT)0 -') accounts 
for the orientational disorder of the ti g spins 2J] ; t^ft m , 
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FIG. 3: (Color online) Orbital polarization p (left) and (right) 
occupied state \0) (| — 9)) for sites 1 and 3 (2 and 4) as a 
function of temperature. Solid line: 300 K (Rn) and 800 K 
(Rl ! 114 ) structures. Dots: orthorhombic structures with half 
(Re) or no (Ro) Jahn- Teller distortion. Pentagons: 2 (full) 
and 4 ( empty) sites CDMFT. Dashes: ideal cubic structure 
(Io). Circles: U — 5 eV. Diamonds: U = 5.5 eV. Triangles: 
U = 6 eV. Squares: U = 7 eV. Crystal-field splitting (meV): 
840 (Ru), 495 (R 6 ), 219 (R ), 168 (Rl° 4 0K ), and (I ). 



in line with experiments We find that even at 

1150 K the system is fully orbitally polarized (p ~ 1). 
On sites 1 and 3, the occupied state is \9) ~ |106°), 
on sites 2 and 4 it is | - 9) ~ | - 106°) (d-type OO); 
\9) is close to the lower crystal- field state obtained from 
LDA (table Q} and in excellent agreement with neutron 
diffraction experiments Q. We find that things hardly 
change when the JT distortion is halved (Kq structure in 
Fig.©. Even for the average 800 K structure (R^) OO 
does not disappear: Although the Jahn- Teller distortion 
is strongly reduced to Sjt = 2.4%, the crystal-field split- 
ting is ~168 meV and the orbital polarization at 1150 K 
is as large as p ~ 0.65, while 9 is now close to 90°. For 
all these structures, orbital order is already determined 
by the distortions via the crystal-field splitting. 

To find the temperature Tkk at which Kugel-Khomskii 
super-exchange drives orbital-order we consider the zero 
crystal-field limit, i.e. the ideal cubic structure, In. The 
e g band- width increases to W e ~ 3.7 eV and for U = 
5 eV the system is a Mott insulator with a tiny gap only 
below T ~ 650 K. We find T KK ~ 650 K, very close 
to the metal- insulator transition (Fig. \S§ . To check how 
strongly Tkk changes when the gap opens, we increase 
U. For U = 5.5 eV we find an insulating solution with 
a small gap of ~ 0.5 eV and Tkk still close to ~ 650 K. 
For U = 6 eV, E g ~ 0.9 eV and T KK ~ 550 K. Even with 
an unrealistically large U = 7 eV, giving E g ~ 1.8 eV, 
Tkk is still as large as ~ 470 K. Thus, despite the small 
gap, Tkk decreases as ~ 1/U, as expected for super- 
exchange. For a realistic U ~ 5 eV, the calculated Tkk ~ 
650 K is surprisingly close to the order-disorder transition 
temperature, Too ~ 750 K, though still much smaller 
than Tjt > 1150 K. The occupied state at site 1 is \9) ~ 



|90°) for all U. 

Such a large Tkk is all the more surprising when com- 
pared with the value obtained for KCuF 3 , Tkk ~ 350 K 
2lT |. For the ideal cubic structure the hopping matrix 
(table © is 4^± z , ~ -t5 m , m >6 m ,3 Z 2-u = & 



m . rn 
i,z±x 



-t 



i,i±y 



—4/4(1 + 2S m x 2_ y 2), and for m ^ m' t, 

Since the effective (after averaging over the di- 
rections of Sj 2 ) hopping integral in LaMnOs, 24/3 ~ 
345 meV is ~ 10% smaller than t ~ 376 meV in KCuF 3 



21j , one may expect a slightly smaller Tkk in LaMnOs, 
opposite to what we find. Our result can, however, be 
understood in super-exchange theory. The KK SE part of 
the Hamiltonian, obtained by second-order perturbation 
theory in t from Eq. ([T]), may be written as 

Hs'e - [ 3T ^ ^ vs (Tin + ?tn 



JsE 



(2) 



where (i, i') x ,y and (i, i'} z indicate near neighboring sites 
along x, y, or z; — (+) refers to the x (y) direction, Tf 
and Tf are pseudospin operators 0, with T z |3z 2 - 1) = 
l/2|3z 2 - 1), T z \x 2 - y 2 ) = -l/2\x 2 - y 2 ). The su- 
perexchange coupling is Jse — (P /U)(w/2), where t 
is the effective hopping integral. In the large U limit 
(negligible J/U and h/U), w - 1 + 4(5? )(5f,) + (1 - 

4{Sf){Sf,))u % 2 l ^/u l 2\, where Sf are the e g spin operators. 
In LaMnOs the e g spins align with the randomly oriented 
t 2g spins, thus t = 24/3, w ~ 2, and J S e ~ 2(2t/3) 2 /U. 
For (4-type order, the classical ground-state is \ ff) ~ |90°), 
in agreement with our DMFT results. In KCUF3, with 
configuration 4^6^, the Hund's rule coupling between e g 
and 42 9 plays no role, i.e. (Sf ) = 0. The hopping inte- 
gral 4 — 4 is indeed slightly larger than in LaMnOs, but 
w ~ 1, a reduction of 50%. Consequently, Jse is reduced 
by ~ 0.6 in KCUF3. For finite J/U and h/U, w is a more 
complicated function, but the conclusions stay the same. 
We verified solving JT]) with LDA+DMFT that also for 
LaMnOs Tkk drops drastically if uf _ a = and h — 0. 

It remains to evaluate the effect of the orthorhombic 
distortion on the transition. For this we perform calcu- 
lations for the system Ro with no Jahn- Teller distortion, 
but keeping the tetragonal and GdFe0 3 -type distortion 
of the 300 K structure. This structure is metallic for 
U = 4 eV; for U = 5 eV it has a gap of — 0.5 eV. We 
find a large polarization already at 1150 K (p ~ 0.45). 
Such polarization is due to the crystal-field splitting of 
about 219 meV, with lower crystal-field states at site 1 
given by |1) ~ \x 2 — y 2 ). Surprisingly the most occupied 
state \9) is close to |1) (6 ~ 180) only at high temper- 
ature (~ 1000 K). The orthorhombic crystal-field thus 
competes with super-exchange, analogous to an external 
field with a component perpendicular to an easy axis. 
On cooling the occupied orbitals rotate to \9) ~ 1 132°) 
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TABLE I: Hopping integrals t*^ , /meV from a site i of type 
1 to a neighboring site i' of type 2 in direction ix + ny + rnz 
for structures Rn, R2°4 K , Ro and Io- The states m,m! are 
|tt) = |x 2 - j/ 2 ) and |0) = |3z 2 - 1). The crystal-field states 
are the eigenvectors of the on-site matrix (I = m = n = 0). 

(see Fig. [3]). This effect of super-exchange, occurs around 
a characteristic temperature T^ K ~ 550 K; still surpris- 
ingly large, but reduced compared to Irk for the ideal 
cubic system Io and much smaller than the experimen- 
tal Tjt > 1150 K. Short-range correlations could reduce 
Trk or rnodify 9. To estimate this effect we perform 
CDMFT calculations; our results (Fig. |3]) remain basi- 
cally unchanged. Thus, electron-phonon coupling is nec- 
essary to explain both the transition temperature and 
the correct occupied orbital \6) ~ |108°). 

In conclusion, we find that T£ K in orthorhombic 
LaMnC>3 is ~ 550 K. We have shown that two ele- 
ments are crucial: the super-exchange mechanism, which 
yields a transition temperature as high as 650 K, and 
the tetragonal plus GdFeC>3-type distortion, which, due 
to the reduced hopping integrals and the competing or- 
thorhombic crystal-field, reduces Tkk to 550 K. Exper- 
imentally, an order-to-disorder transition occurs around 
Too ~ 750 K, but a local Jahn- Teller distortion per- 
sists in the disordered phase up to Tjt > 1150 K. The 
Kugel-Khomskii mechanism alone cannot account for the 
presence of such Jahn- Teller distortions above 550 K 
(T| K <C Tjt)- It also cannot justify the neutron scatter- 
ing estimate 9 = 108°. Thus electron-phonon coupling 
is a crucial ingredient, both for making the Jahn- Teller 
distortions energetically favorable at such high tempera- 
tures and in determining the occupied orbital. 

We acknowledge discussions with D.I. Khomskii and 
A.I. Lichtenstein, and grant JIFF22 on Jugene. 
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